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Abstract 

This  work  presents  a  method  for  extracting  vessel  arrival  times  and  arrival 
processes  from  Automatic  Identification  System  (AIS)  data.  This  work 
employs  the  methodology  presented  by  Mitchell  and  Scully  (2014)  for 
inferring  tidal  elevation  at  the  time  of  vessel  movement  and  calculating  the 
tidal  dependence  (TD)  parameter  to  23  U.S.  port  areas  for  the  years  2012- 
2014.  Tidal  prediction  stations  and  observation  reference  lines  are 
catalogued  for  considered  ports.  AIS  data  obtained  from  the  U.S.  Coast 
Guard,  and  6-minute  tide  predictions,  obtained  from  the  National 
Oceanographic  and  Atmospheric  Administration,  are  used  to  rank  relative 
tidal  dependence  for  arriving  cargo  and  tank  vessel  traffic  in  studied 
locations.  Results  include  relevant  tide  range  and  elevation  threshold 
observations  for  each  year  and  location  studied.  AlS-derived  arrival 
processes,  including  arrival  frequency,  arrival  rate,  and  interarrival  time  are 
visualized  using  several  techniques  with  comparative  discussion  between 
ports  to  highlight  implications  for  understanding  seasonal  traffic  trends  or 
port  resiliency.  The  ports  with  the  highest  and  lowest  TD  value,  Portland, 
ME,  and  Los  Angeles,  CA,  respectively,  are  discussed  with  regard  to  weekly 
arrival  patterns  and  interarrival  time.  Cargo  composition  and  value 
obtained  through  the  Channel  Portfolio  Tool  is  also  considered. 
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1  Introduction 

Background 

The  U.S.  Army  Corps  of  Engineers  (USAGE)  is  responsible  for  planning, 
constructing,  and  maintaining  a  vast  nation-wide  network  of  navigation 
channel  infrastructure  in  coastal  and  riverine  systems.  Environmental 
conditions  vary  widely  across  the  network,  specifically  with  regard  to 
tidally  driven  water  surface  elevations.  Tidal  regimes  include  diurnal, 
semidiurnal,  and  mixed  systems  with  virtually  no  tidally  driven  water 
surface  changes  (e.g.,  on  the  Great  Lakes  and  inland  rivers)  to  tidal  ranges 
that  approach  30  feet  (ft)  (e.g.,  parts  of  Alaska).  Fluctuation  of  the  water 
surface  resulting  from  tide  is  significant  at  projects  that  experience  the 
phenomenon.  Vessel  operators  may  take  advantage  of  tidally  driven  water 
surface  changes  to  sail  with  drafts  larger  than  could  otherwise  be  achieved 
without  additional  tidal  height. 

Recent  examples  of  planning  feasibility  studies  including  harbor  deepening 
at  Savannah,  GA,  Charleston,  SC,  and  the  Port  of  New  York  and  New  Jersey, 
among  others,  demonstrate  that  each  project  considers  the  unique  impacts 
of  tidally  varying  depth  when  planning  for  project  expansion.  Despite  the 
variety  and  importance  of  tidal  fluctuations  within  the  portfolio  of 
navigation  projects,  a  comprehensive  evaluation  of  tidal  influence  across  the 
network  has  not  been  undertaken.  Possible  reasons  for  the  lack  of  a 
comprehensive  evaluation  include  lack  of  access  to  data,  lack  of  an  objective 
and  repeatable  methodology,  and  the  sheer  scope  of  the  problem. 

Objective 

The  objective  of  this  report  is  to  capitalize  on  the  emergence  of  Automatic 
Identification  System  (AIS)  technology  as  a  remote  sensing  tool  for  vessels 
operating  within  the  USAGE  navigation  project  portfolio.  Due  to  the 
expansive  use  of  AIS  by  vessels  within  the  coastal  portion  of  the  navigation 
portfolio,  an  opportunity  exists  to  undertake  an  evaluation  of  navigating 
vessel  behaviors,  including  those  related  to  tidally  driven  water  surface 
fluctuations,  in  a  large  sample  of  USAGE  projects. 

Generally,  AIS  data  cover  most  large  commercial  ships  transiting  coastal 
and  inland  navigation  projects  within  the  USAGE  portfolio.  Goastal  data 
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acquisition  and  dissemination  is  performed  by  the  U.S.  Coast  Guard 
(USCG).  AIS  coverage  for  inland  waterways  is  generally  less  robust  than 
for  coastal  projects.  Inland  data  coverage  is  expanding,  with  collection 
performed  by  the  USAGE.  Once  collected,  data  from  inland  regions  are 
incorporated  into  the  larger  USCG  data  store.  Scully  and  Mitchell  (2015) 
provide  insight  into  potential  uses,  interpretation,  and  availability  of  AIS 
data. 

In  the  course  of  investigating  the  role  of  tide  during  vessel  arrivals,  it  was 
further  recognized  that  due  to  the  nature  of  the  tidal  analysis  performed  in 
this  study,  arrival  processes  can  be  easily  mined  from  AIS  data  as  a 
byproduct  of  the  tidal  analysis  information.  Both  tidal  considerations  and 
traffic  volume  and  frequency  estimation  are  identified  in  USAGE 
engineering  manuals  as  important  navigation  study  inputs  (USAGE  2006). 
Vessel  arrival  processes,  including  arrival  time,  interarrival  time,  and 
arrival  frequency  are  commonly  used  as  inputs  to  USAGE  navigation 
feasibility  study  models.  Arrival  processes  are  normally  derived  from  data 
after  aggregating  a  variety  of  sources,  including  vessel  pilot  and  terminal 
operations  logs,  into  one  dataset.  The  availability  of  a  method  that  relies 
on  AIS  data  allows  for  simplification  and  standardization. 

This  report  is  organized  as  follows.  Chapter  2  describes  the  data  that  was 
used  in  this  investigation  and  how  data  can  be  obtained.  The  methods  for 
performing  tidal  analysis  and  determining  vessel  arrival  processes  are 
described  in  Chapter  3.  Analysis  results  are  presented  in  Chapter  4. 
Chapter  5  contains  discussion  of  the  results  and  their  implications  for 
waterway  managers.  Conclusions  drawn  from  the  investigation  are 
presented  in  Chapter  6. 
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2  Data  Description  and  Sources 

This  investigation  relies  on  data  from  two  primary  sources.  AIS  data  was 
requested  manually  in  full  resolution  from  the  USCG  via  the  agency’s 
Nationwide  AIS  (NAIS)  system.  The  NAIS  system  is  designed  to  enhance 
maritime  domain  awareness  within  continental  United  States  and 
territorial  waters.  Archived  NAIS  data  are  available  from  the  program’s 
historic  data  request  page:  www.navcen.uscg.gov/?pageName=dataReauest& 
dataReauest=aisHistoricalReauestForm.  AIS  data  archives  are  also  available  from 
commercial  vendors  or  may  be  collected  with  the  use  of  an  AIS  receiver 
with  archival  capability.  Water  surface  elevation  information  was 
requested  automatically  through  the  National  Oceanographic  and 
Atmospheric  Administration  (NOAA)  tides  and  currents  applications 
program  interface  (API):  http://tidesandcurrents.noaa.gov/api/. 

AIS  data 

AIS  data  formatting  is  specified  by  the  International  Telecommunication 
Union  (ITU  2014).  The  data  contained  in  position  report  messages  include 
a  time  stamp,  latitude,  longitude,  and  course  over  ground.  Message  type  5 
contains  vessel  particulars,  including  vessel  and  cargo  type  code.  These 
four  data  dimensions  were  used  in  this  study.  The  time  stamp  is  provided 
in  Coordinated  Universal  Time  (UTC).  Position  information  is  provided 
relative  to  the  World  Geodetic  Survey  (WGS)  1984  datum.  Time  stamp  and 
position  information,  both  referred  to  as  dynamic  information,  are 
automatically  generated  by  electronic  navigation  systems  onboard  vessels 
using  AIS  technology.  Position  reports,  contained  in  message  types  1,  2, 
and  3,  are  generated  at  varying  frequencies  that  depend  on  vessel 
behavior,  ranging  from  2  to  180  seconds  (sec)  as  described  in  Table  2.  The 
vessel  and  cargo  type  dimension,  referred  to  as  static  since  it  changes  very 
infrequently,  is  manually  generated.  This  is  a  two-digit  code,  described  in 
the  ITU  technical  specification.  Several  authors,  including  the  commonly 
cited  work  by  Harati-Mokhtari  et  al.  (2007),  describe  considerations  that 
must  be  given  to  ensure  AIS  data  quality  is  sufficient  for  specific  use  cases. 

Water  surface  elevation  data 

NOAA  predicts  water  surface  elevation  data  at  several  frequencies  for 
locations  that  are  influenced  by  the  tide.  NOAA  also  collects  water  surface 
elevation  information,  which  can  be  compared  to  predictions.  This  study 
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relied  on  water  surface  predictions  generated  at  6 -minute  (min)  intervals 
for  major  stations  within  the  ports  of  interest.  Water  surface  elevation  data 
is  referenced  to  mean  lower  low  water  (MLLW)  elevation  datum,  and 
prediction  times  are  in  UTC.  Water  surface  predictions,  instead  of  verified 
water  surface  elevations,  are  used  for  three  general  reasons.  First,  it  is 
most  likely  that  tide-reliant  vessels  planning  transits  will  rely  on 
predictions  to  plan  transits,  and  verified  information  will  be  unavailable 
for  planning.  Second,  verified  information  includes  the  effects  of 
atmospheric  and  hydrologic  influence  that  are  most  likely  well  beyond  the 
ability  of  NOAA  to  predict  accurately.  Third,  operational  limitations 
including  loss  of  tidal  station  measurement  function  may  result  in  data 
gaps  within  the  study  period. 
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3  Methods 

Scully  and  Mitchell  (2013)  and  Mitchell  and  Scully  (2014)  outline  a 
method  for  interpretation  of  water  vessel  surface  elevation  information  at 
the  time  of  vessel  arrival.  This  technical  report  expands  the  level  of 
detailed  methodology  provided  by  the  authors. 

Observation  location  selection 

A  reference  observation  line  is  required  to  investigate  tidal  influence.  A 
line  crossing  a  waterway,  defined  by  two  latitude  and  longitude  pairs, 
should  be  selected  at  the  location  of  interest.  This  line  will  be  used  to  count 
the  number  of  vessels  crossing  the  line  and  to  document  the  time  of  vessel 
crossing.  Both  are  key  methodological  inputs. 

Observation  reference  lines  should  be  carefully  selected  and  based  on 
known  navigation  channel  locations.  Alternatively,  reference  location 
selection  maybe  made  using  geographical  features  such  as  coastal  inlets  or 
exploratory  plotting  of  AIS  data  to  identify  vessel  position  report  clusters. 
AIS  information  is  carried  over  very  high  frequency  (VHF)  radio  signals, 
which  are  limited  to  line-of-sight  transmission.  Geographical  features, 
including  bridges,  valleys,  dense  urban  structures,  or  anything  conceivably 
capable  of  interfering  with  AIS  transmission  or  reception,  should  be  avoided 
when  selecting  reference  locations.  Sample  reference  lines  and  tidal  stations 
used  in  this  study  are  shown  for  Seattle  and  Tacoma,  WA;  New  York  and 
New  Jersey;  Los  Angeles  and  Long  Beach,  CA;  and  Corpus  Christi,  TX,  in 
Figure  1.  A  reference  line  for  each  port  is  shown  in  Appendix  A. 

Consideration  should  be  given  to  the  types  of  vessels  that  are  likely  to  use 
the  waterway  being  investigated.  Reference  lines  should  be  selected  to 
capture  the  vessel  population  of  interest.  Challenges  may  be  encountered 
in  locations  with  complex  channel  network  components  such  as  loops  or 
branches.  Locations  with  clustered  terminals  and  short-range  shipping 
routes  may  also  require  further  assessment  of  reference  selection. 
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Figure  1.  Example  tide  station  and  reference  line  arrangement. 
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This  study  prioritized  capture  of  the  greatest  fraction  possible  of  large 
commercial  vessel  traffic  by  focusing  on  coastal  inlets  to  port  areas.  The 
observation  reference  lines  used  in  this  study  are  listed  in  Table  i. 


Table  1.  Port  areas,  reference  line  coordinates,  and  NOAA  (2013a)  tide  stations. 


Port  Area  Entrance 

Tide  Station 
Number 

Line  Beginning 
Coordinate 

Line  Ending  Coordinate 

Anacortes,  WA 

9444900 

(-122.7080,  48.4792) 

(-122.7109,  48.5482) 

Baltimore,  MD 

8574680 

(-76.5651,  39.2375) 

(-76.5493,  39.2562) 

Boston,  MA 

8443970 

(-71.0064,  42.3192) 

(-71.0080,  42.3504) 

Charleston,  SC 

8665530 

(-79.8546,  32.7592) 

(-79.8703,  32.7363) 

Columbia  River,  OR 

9439040 

(-124.0860,  46.2644) 

(-124.0750,  46.2331) 

Delaware  Bay,  DE 

8537121 

(-75.4122,  39.2619) 

(-75.3200,  39.3253) 

Freeport,  TX 

8772447 

(-95.2930,  28.9296) 

(-95.2898,  28.9329) 

Honolulu,  HI 

1612340 

(-157.8656,  21.2955) 

(-157.8741,  21.3020) 

Jacksonville,  FL 

8720218 

(-81.4214,  30.4094) 

(-81.4216,  30.3976) 

Long  Beach,  CA 

9410660 

(-118.1932,  33.7230) 

(-118.1726,  33.7226) 

Los  Angeles,  CA 

9410660 

(-118.2544,  33.7080) 

(-118.2370,  33.7134) 
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Port  Area  Entrance 

Tide  Station 
Number 

Line  Beginning 
Coordinate 

Line  Ending  Coordinate 

Mobile,  AL 

8737048 

(-88.0447,  30.6889) 

(-88.0301,  30.6923) 

New  Haven,  CT 

8465705 

(-72.9348,  41.2603) 

(-72.9018,  41.2480) 

New  York  and  New  Jersey 

8531680 

(-74.0554,  40.6036) 

(-74.0348,  40.6096) 

Norfolk,  VA 

8638610 

(-76.3150,  36.9513) 

(-76.3595,  37.0049) 

Pascagoula,  MS 

8741533 

(-88.5698,  30.3296) 

(-88.5033,  30.3202) 

Port  Everglades,  FL 

8723214 

(-80.1050,  26.0960) 

(-80.1055,  26.0925) 

Portland,  ME 

8418150 

(-70.2390,  43.6662) 

(-70.2238,  43.6521) 

San  Francisco  Bay,  CA 

9414290 

(-122.4770,  37.8094) 

(-122.4800,  37.8268) 

San  Juan,  PR 

9755371 

(-66.1361,  18.4718) 

(-66.1234,  18.4708) 

Savannah,  GA 

8670870 

(-80.8863,  32.0150) 

(-80.8913,  32.0684) 

Seattle-Tacoma,  WA 

9447130 

(-122.4064,  47.6535) 

(-122.5190,  47.6602) 

Tampa,  FL 

8726384 

(-82.7348,  27.6224) 

(-82.7394,  27.5337) 

Interpret  water  surface  threshold  elevations 

This  study  used  high-,  mid-,  and  low-tide  segments  of  the  tidal  cycle  (Scully 
and  Mitchell  2013).  High  tide  and  low  tide  are  defined  as  the  upper  and 
lower  time-based  quartile  of  water  surface  elevations,  respectively.  Mid-tide 
is  defined  as  the  time-based  interquartile  range  of  tidal  elevations.  Each 
year’s  record  of  predicted  tidal  elevation  (87,600  predictions  at  6  min 
intervals  for  a  standard  year)  was  ordered  by  elevation.  The  25th  quartile 
elevation  is  the  threshold  between  low  and  mid  tide  and  represents  the 
elevation  below  which  tidal  elevations  occurred  25%  of  the  year.  Similarly, 
the  75th  quartile  elevation  is  the  threshold  between  mid  and  high  tide  and 
represents  the  elevation  above  which  tidal  elevations  occurred  25%  of  the 
year.  The  mid  tides  occurred  50%  of  the  year  between  the  high  and  low 
thresholds.  Thresholds  for  each  location  are  determined  independently  for 
each  of  the  years  2012,  2013,  and  2014. 

Vessel  filtering,  transit  generation,  and  water  surface  at  time  of 
crossing 

AIS  data  includes  many  dimensions  for  data  filtering  (ITU  2014;  Scully 
and  Mitchell  2015).  Filtering  for  this  study  was  applied  to  static  and 
dynamic  data  components.  Data  were  initially  received  from  the  USCG  as 
a  collection  of  i-month-increment,  comma-separated  value  files,  organized 
by  location.  The  data  were  processed  using  the  Python  programming 
language  (van  Rossum  and  Drake  2001),  the  Pandas  data  analysis  package 
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(McKinney  2012),  or  the  developmental  version  of  USAGE  AIS  Analysis 
Package  (AISIAP)  software. 

Individual  vessels  are  defined  in  this  study  as  those  having  a  unique 
Maritime  Mobile  Service  Identifier  (MMSI)  number.  Several  authors  have 
written  about  the  potential  complications  of  this  approach  (Harati- 
Mokhtari  et  al.  2007).  By  eliminating  or  verifying  duplicate,  vague,  or 
incorrect  MMSI  numbers,  most  of  these  complications  can  be  mitigated. 
For  high-accuracy  applications,  comparing  MMSI  data  to  verified 
authoritative  vessel  information  is  recommended. 

Vessels  were  filtered  in  this  study  based  on  the  static  “ship  and  cargo  type” 
dimension,  contained  in  Message  5  (again,  subject  to  human  error 
implications).  Vessels  with  ship-  and  cargo-type  codes  that  began  with 
only  “7”  (cargo  ships)  or  “8”  (tankers)  were  included.  Vessels  were  also 
filtered  on  the  dynamic  temporal,  spatial,  and  heading  dimensions, 
primarily  derived  from  Messages  1,  2,  or  3. 

Vessels  using  Class  A  AIS  transceivers  report  their  position  according  to 
the  reporting  frequencies  defined  by  the  ITU  (2014)  and  shown  in  Table  2. 
A  simple  approach  for  identifying  unique  transits  is  to  identify  gaps  in  the 
position  report  sequence  that  exceed  normal  report  intervals. 


Table  2.  Class  A  AIS  vessel  reporting  Intervals  (ITU  2014,  Table  1). 


Ship's  Dynamic  Conditions 

Nominal  Reporting  Interval 

Ship  at  anchor  or  moored  and  not  moving  faster  than  3  knots 

3  min 

Ship  at  anchor  or  moored  and  moving  faster  than  3  knots 

10  sec 

Ship  0-14  knots 

10  sec 

Ship  0-14  knots  and  changing  course 

3  1/3  sec 

Ship  14-23  knots 

6  sec 

Ship  14-23  knots  and  changing  course 

2  sec 

Ship  >  23  knots 

2  sec 

Ship  >  23  knots  and  changing  course 

2  sec 

Transits  were  defined  by  identifying  gaps  in  the  position  report  sequence 
that  exceed  360  sec  (Scully  and  Mitchell  2013).  This  duration  is  chosen  to 
ensure  that  vessels  reporting  their  position  at  the  lowest  specified  frequency 
(i.e.,  a  “ship  at  anchor  or  moored  and  not  moving  faster  than  3  knots”)  will 
have  at  least  one  position  identified  if  it  is  reporting  normally  and  its  signal 
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is  unobstructed.  These  vessels  are  not  generally  of  interest  when  assessing 
water  surface  elevation  during  transit.  Ships  moving  with  speed  between  o 
and  14  knots  report  at  the  next  most  frequent  interval,  one  report  every  10 
sec.  Vessels  moving  with  speeds  in  this  range  will  have  ample  position 
reports  to  construct  a  transit.  Vessel  reports  that  have  been  smoothed  to  a 
frequency  lower  than  one  report  per  3  min  will  require  a  longer  gap  length 
or  different  transit  generation  method  to  reliably  treat  the  data. 

Vessel  transits  are  limited  by  spatial  filters.  First,  a  5,000  ft  buffer  was 
applied  to  either  side  of  each  observation  reference  line.  Minimum 
bounding  envelope  geometry  was  applied  to  buffered  areas  to  create 
regularly  arranged  rectangles  for  filtering.  Vessels  crossing  each  reference 
line  are  expected  to  transit  on  the  order  of  10,000  ft  across  the  local  study 
area. 

The  shortest  report  duration  is  specified  as  one  report  every  2  sec  for 
vessels  transiting  at  speeds  greater  than  23  knots.  A  vessel  moving  at 
23  knots  (38.8196  ft/sec)  will  sail  10,000  ft  in  approximately  258  sec  and 
report  its  position  approximately  129  times.  This  buffer  distance  ensures 
that  transiting  vessels  among  the  population  of  interest  will  generate  at 
least  two  position  reports  within  each  study  area,  ensuring  transit 
generation. 

Vessels  were  classified  as  either  transiting  inbound  or  outbound  based  on 
the  heading  data  dimension.  To  facilitate  this  classification,  an  inbound 
direction  vector  was  arbitrarily  defined  as  being  approximately  normal  to 
each  observation  reference  line  and  pointing  away  from  the  open  ocean 
and  toward  the  port  area  of  interest.  Vessels  were  classified  as  inbound  if 
the  course  over  ground  values  of  each  transit  were  within  +/-  90°  of  the 
cartographic  direction  of  the  inbound  definition  vector,  and  outbound 
otherwise.  For  example,  if  the  inbound  direction  vector  is  defined  to  have  a 
cartographic  heading  of  300°,  then  inbound  vessels  are  those  with  course 
over  ground  ranging  from  210°  to  360°  and  from  0°  to  30°.  Only  inbound 
vessels  were  included  in  this  study. 

To  summarize  the  filtering  scheme,  only  unique  tanker  and  cargo  vessels 
(based  on  unique  identifying  information  and  ship  and  cargo  type  code 
contained  in  AIS),  heading  away  from  the  open  ocean  and  toward  selected 
ports  (based  on  AIS  embedded  course  over  ground  information)  were 
considered.  Vessels  were  selected  from  conservative  buffer  regions  around 
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observation  reference  lines  chosen  to  capture  deep-draft  vessel  traffic. 
Analysis  was  performed  for  each  location  discretely  in  2012,  2013,  and 
2014.  Aside  from  this  filtering,  data  were  not  thoroughly  controlled  for 
quality  owing  to  the  large  quantity  of  data.  This  induces  the  risk  that  a 
fraction  of  the  vessel  population  is  missing  or  incorrectly  classified  within 
the  AIS  record.  In  all,  over  120,000  vessel  transits  were  generated  to 
assess  tidal  influence. 

Linear  interpolation  is  used  to  determine  the  time  at  which  vessels  cross 
observation  reference  lines.  The  nearest  position  report  to  either  side  of 
the  reference  line  is  selected,  and  the  time  difference  between  reports  is 
calculated.  The  distance  between  these  two  positions,  and  the  distance 
from  the  earlier  position  report  to  the  observation  reference  line,  are  also 
determined.  The  vessel-crossing  offset  is  computed  as  the  proportion  of 
the  reference  line  distance  to  the  total  distance.  The  offset  is  then 
multiplied  by  the  time  difference  between  position  reports  and  added  to 
the  time  of  the  earlier  position  report,  yielding  the  time  of  vessel  crossing. 

Interpret  water  surface  elevation  at  time  of  vessel  crossing 

The  water  level  at  the  time  each  vessel  crossed  a  reference  line  is  linearly 
interpolated  from  the  water  level  record  using  the  time  of  vessel  crossing 
and  assigned  to  that  transit.  Vessels  in  transit  generally  move  faster 
laterally  than  the  tide  can  rise  or  fall  over  short  periods  of  time.  The  most 
extreme  tides  in  the  world,  at  the  Bay  of  Fundy,  have  been  documented  to 
experience  changes  in  water  surface  elevation  of  56  ft  in  approximately 
12.42  hours  (hr)  (NOAA  2013b).  The  water  surface  elevation  changes  at  a 
rate  of  1.3  x  i0‘3  ft/sec  (7.4  x  io-4  knots).  A  vessel  reporting  its  position 
every  3  min  would  observe  a  change  in  water  surface  elevation  of 
approximately  0.3  ft  at  that  location.  Since  most  vessels  of  interest  to  this 
study  report  their  position  every  10  sec,  and  the  tidal  variations  at  study 
locations  are  much  less  than  the  example  (Seattle  had  the  greatest  tidal 
range  at  16.71  ft),  it  is  assumed  the  resulting  errors  in  water  surface 
elevation  at  time  of  vessel  crossing  are  negligible. 

It  was  further  assumed  that  errors  in  water  surface  interpolation  resulting 
from  hydraulic  friction  loss  over  the  distance  between  prediction  location 
and  observation  location  were  negligible.  For  instance,  the  tidal  prediction 
station  at  the  port  of  Anacortes,  WA,  provides  high-  and  low-tide 
predictions,  based  on  the  harmonic  tidal  observing  station  9444900  at 
Port  Townshend,  WA,  but  does  not  provide  6  min  water  level  predictions. 
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The  Port  Townshend  station  is  approximately  28  miles  from  the 
observation  line  used  to  analyze  Anacortes  traffic.  The  Anacortes 
prediction  station  provides  a  time  and  elevation  offset  to  the  harmonic 
station.  High  tides  at  the  Anacortes  are  predicted  to  occur  22  min  later 
than  at  Port  Townshend,  and  the  tidal  height  is  approximately  96%  of  the 
high-tide  height.  Low  tide  occurs  approximately  33  min  later  at  Anacortes 
than  Port  Townshend,  and  the  low-tide  heights  are  approximately  equal. 

The  difference  in  high-tide  height  between  Anacortes  and  Port  Townshend 
is  less  than  6  inches  (in.),  or  3.9%  of  the  observed  tidal  range.  The  diurnal 
tide  range  for  Port  Townshend  is  9  ft,  and  the  tidal  period  is  approximately 
12.42  hr,  meaning  a  normal  tide  celerity  of  0.72  ft  per  hr.  A  33  min  lag 
would  result  in  approximately  5  in.  in  tide  differential.  The  combined  tide 
height  and  time  lag  errors  amount  to  approximately  1  ft  of  the  observed 
average  tide  range  of  12.98  ft,  which  amounts  to  approximately  7.7%  error. 
Because  the  observation  reference  line  for  Anacortes  is  closer  to  the  Port 
Townshend  harmonic  station  than  the  Anacortes  prediction  station,  it  can 
be  expected  that  induced  errors  will  be  smaller. 

Calculate  tide  sector  traffic  percentages  and  tidal  dependence  metric 

After  assigning  the  water  level  to  each  observed  transit,  compare  water 
level  transits  to  tide  thresholds  to  apportion  traffic  to  tidal  cycle  segments. 
The  proportion  of  vessels  assigned  water  surface  elevations  below  the  low- 
tide  threshold  is  T25.  The  proportion  of  vessels  transiting  above  the  high- 
tide  threshold  is  T75.  T50,  the  portion  of  traffic  transiting  during  mid-tide, 
can  be  calculated  as  1  -  (T25  +T75). 

The  tidal  dependence  (TD)  metric  (Mitchell  and  Scully  2014)  is  calculated 
as 


TD={T„-T^)/T„  (1) 

The  refined  tidal  dependence  metric  is  intuitive.  It  is  formulated  such  that 
vessel  traffic  uniformly  distributed  across  the  tidal  cycle  will  result  in  a  TD 
value  of  unity.  Traffic  regimes  with  dominant  high-tide  traffic  will  result  in 
positive  TD  values  whereas  low-tide  dominated  ports  will  be  negative.  Mid¬ 
tide  traffic  acts  to  scale  the  value— as  the  mid-tide  proportion  of  transits 
increases  the  TD  parameter  approaches  zero.  As  the  mid-tide  proportion  of 
transits  decreases,  the  balance  of  high-tide  vs.  low-tide  traffic  becomes  more 
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apparent,  driving  the  value  toward  +/-  oo  asymptotically  until  the  mid-tide 
portion  equals  o. 

Arrival  process  mining 

Linear  interpolation  of  the  time  of  vessel  arrival  was  demonstrated  for 
single  vessels  in  the  previous  section  “Vessel  filtering,  transit  generation, 
and  water  surface  at  time  of  crossing.”  The  time  of  vessel  arrival  is  a 
critical  input  to  port  feasibility  studies  and  is  a  fundamental  input  to 
arrival  rate  and  interarrival  process  mining.  Vessel  arrival  time  is  a 
required  input  to  models  that  seek  to  quantify  the  impacts  to  changes  in 
navigation  infrastructure.  Aggregating  individual  arrival  times  for  a  vessel 
population  at  the  reach  or  port  level  results  in  the  arrival  frequency 
distribution  of  the  navigation  feature,  which  describes  larger  operations. 

Arrival  processes  of  general  interest  to  USAGE  navigation  planners  and 
operators  alike  include  the  arrival  rate  and  the  interarrival  time  of  vessels 
calling  in  a  port  or  reach.  The  arrival  rate  is  a  measure  of  vessel  arrivals 
per  unit  time.  This  measure  is  derived  by  dividing  the  number  of  arrivals 
by  the  unit  time  of  interest.  Vessels  per  day  or  month  are  commonly  used 
as  benchmark  indicators  of  vessel  activity.  Long-term  averages  of  vessel 
arrival  rates  can  be  compared  to  short-term  arrival  rates  to  identify  peaks 
and  lulls  in  vessel  activity.  Interarrival  time  is  calculated  by  ordering  vessel 
arrivals  chronologically  and  calculating  the  time  between  arrivals  to  find 
the  distribution. 


ERDC/CHLTR-17-2 


13 


4  Results 

The  methods  described  in  the  preceding  section  facilitate  tidal  analysis 
and  arrival  process  extraction  between  geographically  separated  locations 
with  differing  tidal  patterns.  Table  3  summarizes  the  overall  description  of 
investigated  port  areas  and  includes  the  distance  from  tide  stations  to 
observation  lines  and  observed  tidal  patterns.  Values  averaged  over  the 
years  2012,  2013,  and  2014  include  mean  threshold  elevations,  tidal  range, 
and  the  number  of  arriving  vessels  included  in  the  analysis. 


Table  3.  Port  area  tidal  information. 


Port  Area  Entrance 

Distance 
from  Tide 
Station 
(miles) 

Tidal  Pattern 

Mean  25th 
Percentile 
Elevation  (ft) 

Mean  75th 
Percentile 
Elevation  (ft) 

Avg.  Tide 

Range 

(ft) 

Avg. 

Number. 

Arriving 

Vessels 

Anacortes,  WA 

28.2 

mixed 

3.19 

7.05 

12.98 

473 

Baltimore,  MD 

1.8 

semi-diurnal 

0.43 

1.16 

2.71 

953 

Boston,  MA 

2.8 

semi-diurnal 

1.99 

8.41 

14.43 

719 

Charleston,  SC 

4.3 

semi-diurnal 

1.11 

4.69 

8.29 

1978 

Columbia  River,  OR 

15.3 

mixed 

2.43 

6.73 

12.53 

1572 

Delaware  Bay,  DE 

2.8 

semi-diurnal 

1.07 

4.92 

8.86 

2161 

Freeport,  TX 

1.1 

mixed 

0.60 

1.37 

3.12 

797 

Honolulu,  HI 

0.7 

mixed 

0.36 

1.25 

3.11 

1014 

Jacksonville,  FL 

0.7 

semi-diurnal 

0.94 

4.02 

7.50 

1541 

Long  Beach,  CA 

5.1 

mixed 

1.69 

3.99 

8.88 

2003 

Los  Angeles,  CA 

1.6 

mixed 

1.69 

3.99 

8.88 

2320 

Mobile,  AL 

1.3 

diurnal 

0.43 

1.20 

2.98 

690 

New  Haven,  CT 

2.1 

semi-diurnal 

1.26 

5.37 

9.13 

138 

New  York  and  New  Jersey 

9.8 

semi-diurnal 

1.03 

4.11 

7.78 

4608 

Norfolk,  VA 

2.3 

semi-diurnal 

0.52 

2.17 

4.11 

3989 

Pascagoula,  MS 

3.3 

diurnal 

0.36 

1.14 

3.00 

764 

Port  Everglades,  FL 

25.1 

semi-diurnal 

0.44 

1.76 

3.66 

3083 

Portland,  ME 

0.7 

semi-diurnal 

1.75 

8.11 

13.84 

269 

San  Francisco  Bay,  CA 

1.1 

mixed 

1.87 

4.48 

8.94 

3143 

San  Juan,  PR 

1.2 

mixed 

0.41 

1.12 

2.62 

1233 

Savannah,  GA 

0.8 

semi-diurnal 

1.49 

6.12 

10.48 

2471 

Seattle-Tacoma,  WA 

6.8 

mixed 

4.30 

9.36 

16.71 

1821 

Tampa,  FL 

11.6 

mixed 

0.71 

1.64 

3.75 

1007 
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Tidal  analysis 

Figure  2  shows  the  average  of  2012,  2013,  and  2014  tidal  segment  traffic 
percentages  for  each  studied  port,  ordered  by  increasing  low-tide  traffic 
percentage.  Portland,  ME,  had  the  greatest  percentage  of  vessel  traffic, 
33.5%,  moving  above  the  high-tide  threshold.  New  Haven,  CT,  and 
Anacortes,  WA,  stand  out  as  having  a  relatively  high  portion  of  high- tide 
traffic.  Los  Angeles,  CA,  had  the  lowest  percentage  of  vessel  traffic,  20.1%, 
transiting  above  the  high-tide  threshold.  Los  Angeles,  CA,  and  Long  Beach, 
CA,  both  stand  out  as  having  a  relatively  high  proportion  (29%)  of  traffic 
arriving  below  the  low-tide  threshold.  Portland  had  the  lowest  fraction  of 
traffic  arriving  below  the  low-tide  threshold,  at  13%. 


Figure  2.  Average  2012-2014  low-,  high-,  and  mid-tide  traffic  distributions. 
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Average  traffic  percentages  and  TD  values  are  summarized  for  each  port  in 
Table  4.  The  fraction  of  traffic  moving  on  the  ends  of  the  tidal  range  for 
Portland,  ME,  is  46%,  compared  to  Los  Angeles,  which  is  49%.  While  these 
two  ports  show  the  greatest  tidal  imbalances,  they  are  not  the  ports  with 
the  greatest  fraction  of  traffic  moving  above  or  below  the  respective  high- 
and  low-tide  thresholds.  Anacortes  saw  the  highest  proportion  of  traffic,  at 
54%,  move  outside  of  mid-tide.  Boston  saw  the  lowest  fraction  of  traffic 
moving  outside  of  mid-tide,  at  43%.  The  variable  range  of  traffic  in  point 
terms  was  similar  for  all  segments:  16  for  low  tide,  11  for  mid  tide,  and  13 
for  high  tide. 
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Table  4.  Average  2012-2014  T25,  T50,  T75,  and  TD  values. 


Port  Area  Entrance 

T25  (mean) 

T50  (mean) 

T75  (mean) 

TD  (mean) 

Anacortes,  WA 

0.22 

0.46 

0.32 

0.24 

Baltimore,  MD 

0.25 

0.48 

0.26 

0.02 

Boston,  MA 

0.16 

0.57 

0.28 

0.21 

Charleston,  SC 

0.25 

0.49 

0.26 

0.01 

Columbia  River,  OR 

0.25 

0.53 

0.23 

-0.01 

Delaware  Bay,  DE 

0.23 

0.50 

0.27 

0.10 

Freeport,  TX 

0.23 

0.50 

0.27 

0.07 

Honolulu,  HI 

0.24 

0.50 

0.26 

0.05 

Jacksonville,  FL 

0.22 

0.52 

0.26 

0.07 

Long  Beach,  CA 

0.29 

0.50 

0.21 

-0.17 

Los  Angeles,  CA 

0.29 

0.51 

0.20 

-0.18 

Mobile,  AL 

0.23 

0.50 
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Figure  3  shows  the  TD  value  of  each  port  in  descending  order.  This  plot 
can  be  interpreted  to  identify  where  preference  for  traffic  moving  in  a 
particular  tide  segment  exists.  Portland  traffic  can  be  said  to  demonstrate 
a  preference  for  high  tide  whereas  Los  Angeles  can  be  said  to  exhibit  low- 
tide  preference.  Charleston,  SC,  and  Columbia  River,  OR,  can  be  said  to 
show  weak  tidal  preference.  It  is  interesting  to  note  that  Boston  and 
Anacortes,  the  ports  with  the  least  and  most  traffic  moving  in  any  tidal 
segment,  respectively,  both  show  a  relatively  high  preference  for  high-tide 
transits. 
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Figure  3.  Average  (2012-2014)  TD  value  for  inbound  cargo  and  tanker  vessels. 
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Arrival  processes 

Tlie  number  of  vessels  arriving  at  each  port  per  day,  calculated  as  the  total 
number  of  observed  arrivals  per  year  divided  by  the  number  of  days  in  that 
year,  provides  an  indication  of  potential  traffic  congestion  within  the  port. 
Arrival  process  mining  from  AIS  data  results  in  the  arrival  rates  displayed 
in  Figure  4.  The  majority  of  studied  ports  have  fewer  than  10  vessel 
arrivals  per  day.  New  York  and  New  Jersey  had  the  highest  arrival  rate 
with  an  average  of  12.6  vessels  per  day.  New  Haven  had  the  lowest  arrival 
rate  (0.4  arrivals  per  day).  The  ports  of  Los  Angeles  and  Long  Beach  had  a 
combined  arrival  rate  of  11.9  vessels  per  day,  ranking  it  as  the  second 
busiest  port  area  in  the  study,  ahead  of  Norfolk  with  a  rate  of  10.9  arrivals 
per  day. 
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Figure  4.  Daily  arrivals  at  studied  ports,  2012-2014. 

Average  Arrivals  Per  Day,  2012-2014 


Tampa,  FL 
Seattle  &  Tacoma,  WA 
Savannah,  GA 
San  Juan,  PR 
San  Francisco,  CA 
Port  Everglades,  FL 
Portland,  ME 
Pascagoula,  MS 
New  York  &  New  Jersey 
Norfolk,  VA 
New  Haven,  CT 
Mobile,  AL 
Los  Angeles,  CA 
Long  Beach,  CA 
Jacksonville,  FL 
Honolulu,  HI 
Freeport,  TX 
Delaware  Bay 
Columbia  River,  OR 
Charleston,  SC 
Boston,  MA 
Baltimore,  MD 
Anacortes,  WA 


5.0 


■  8.6 
8.4 


■  0.4 

1.9 

2.0 

1.3 

10.9 


0.0 


2.0 


4.0 


6.0  8.0 
Arrivals  per  Day 


10.0 


12.0 


14.0 


Figure  5  shows  the  relative  percentage  of  vessel  arrivals  per  day  of  the  week 
at  Los  Angeles,  CA,  and  Portland,  ME,  the  two  most  tidally  influenced  ports 
identified  through  this  analysis.  The  frequency  distributions  are  similar, 
generally  declining  from  Monday  through  Sunday.  This  observation  was 
common  in  many  ports  and  suggests  that  the  previous  arrivals-per-day 
measure  under-reports  arrivals  Monday  through  Friday  and  over-reports 
arrivals  on  Saturday  and  Sunday.  Los  Angeles  traffic  was  observed  to  peak 
on  Monday  and  Wednesday  with  17%  of  traffic  arriving  each  of  those  days. 
Portland  traffic  was  observed  to  arrive  with  peak  frequencies  on  Monday, 
Thursday,  and  Friday,  each  with  approximately  15%  of  traffic  on  each  of 
those  days.  The  least  frequent  arrival  day  in  both  cases  was  Sunday,  with 
arrivals  of  10%  and  13%  of  traffic  arriving  that  day  in  Los  Angeles  and 
Portland,  respectively. 
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Figure  5.  Vessel  arrival  frequencies  are  similarly  distributed  at  Los  Angeles,  CA,  and  Portland,  ME. 
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Greater  detail  of  vessel  arrival  patterns  is  available  from  the  AlS-derived 
vessel  arrival  data.  Figure  6  shows  the  number  of  vessels  observed  to 
arrive  at  the  ports  of  Los  Angeles,  CA,  and  Portland,  ME,  in  2012  by  day  of 
the  week  and  hour  of  the  day.  It  is  evident  from  this  mapping  that 
seemingly  similar  traffic  distributions  have  very  different  daily  patterns. 
Portland  demonstrates  weak  clustering  of  vessel  arrivals  whereas  Los 
Angeles  demonstrates  strong  arrival  clusters  Monday  through  Saturday, 
centered  around  1200  (UTC)  and  2200  (UTC). 


Figure  7  shows  the  interarrival  times  for  cargo  and  tanker  vessels  calling  at 
Los  Angeles,  CA,  and  Portland,  ME,  binned  in  1  hr  increments.  The  time 
between  arrivals  at  Los  Angeles  follows  a  negative  exponential  distribu¬ 
tion.  Fully  35%  of  arrivals  at  Los  Angeles  are  followed  by  another  arrival 
within  1  hr.  These  arrival  processes  can  be  informative  when  drawing 
conclusions  or  making  decisions  related  to  vessel  traffic  within  a  port.  For 
instance,  interarrival  time  distributions  may  be  used  to  inform  discussions 
regarding  the  time  required  to  recover  from  traffic  disruptions. 
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Figure  6.  Daily  vessel  arrivals  by  hour  for  Portland  and  Los  Angeles,  2012. 
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Figure  7.  Interarrival  times  show  differing  levels  of  congestion  at  Los  Angeles,  CA,  and  Portland,  ME,  in  2012. 
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Monthly  arrival  rates  can  be  mined  from  AIS  data  to  help  identify  larger 
traffic  patterns  if  they  exist  within  a  port.  For  example,  this  analysis  of  AIS 
data  suggests  a  spring  peak  in  vessel  arrivals  in  Tampa  in  2012-2014, 
shown  in  Figure  8(a),  and  a  summer  peak  in  arrivals  in  the  Columbia  River 
in  2012  and  2013,  shown  in  Figure  8(b).  The  frequency  range  of  monthly 
arrivals  may  serve  as  a  metric  for  variability  of  port  traffic.  The  Columbia 
River  had  the  greatest  range  of  arrival  frequency  (13.3%),  while  Seattle  and 
Tacoma,  shown  in  Figure  8(c),  had  the  narrowest  range  (2.0%). 

A  complete  set  of  results  for  each  port  is  provided  in  Appendix  A.  Results 
for  each  port  include 

•  general  arrangement  of  the  tide  station  and  observation  reference  line 

•  tidal  station  information 

•  yearly  summary  for  water  surface  thresholds  and  observed  tidal  range 

•  yearly  summary  for  number  of  vessel  arrivals,  vessel  interarrival  time, 
and  arrival  frequency 

•  yearly  summary  for  tide  segment  traffic  percentages 

•  yearly  distributions  of  vessel  arrival  by  water  surface  elevation. 
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Figure  S.Monthly  arrivals  at  (a)  Tampa  and  (b)  Columbia  River  showing  seasonal  increases  in 
vessel  arrivals.  Columbia  River  had  the  widest  range  of  monthly  arrival  frequencies,  (c)  Seattle 
and  Tacoma  had  the  narrowest  range  of  arrival  frequencies 
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(c) 
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5  Discussion 

Prior  to  the  availability  of  AIS  data,  a  multiport  analysis  of  this  scope  and 
detail  was  impractical.  USAGE  typically  analyzes  vessel  movement  data 
provided  from  third  parties  such  as  shipping  terminal  operators  or  harbor 
pilots  in  the  course  of  feasibility  studies  when  considering  harbor 
improvement.  These  sources  may  report  vessel  movement  only  as  they 
apply  to  their  operations.  For  example,  a  terminal  operator  may  report  the 
time  a  vessel  arrives  but  would  have  limited  information  regarding  the  path 
taken  by  a  vessel  or  the  duration  of  its  transit.  Similarly,  pilots  may  log  only 
the  time  they  embark  a  vessel.  Data  generally  do  not  refer  to  the  physical 
location  of  a  vessel  in  transit  but  instead  refer  to  the  times  of  arrival  or 
departure  at  a  particular  landmark.  Each  port  has  many  potential  sources, 
which  must  be  consulted  to  obtain  a  complete  picture  of  harbor  operations. 
There  is  no  standard  format  for  compiling  or  reporting  this  information, 
which  must  be  painstakingly  gathered  and  organized  for  detailed  analysis. 

Recently,  the  Channel  Portfolio  Tool  (CPT)  (Mitchell  2012)  has  improved 
the  ability  to  investigate  the  use  of  navigation  channels  in  support  of 
commodity  movement.  This  tool  relies  on  proprietary  dock-level 
waterborne  commerce  data  and  uses  spatial-join  and  shortest-path 
algorithms  to  aggregate  and  attribute  cargo  tonnage  and  dollar-value 
totals  to  the  respective  transited  reaches.  However,  since  the  tool  is  based 
on  annualized  reported  statistics,  it  does  not  incorporate  the  actual 
behavior  of  vessels  while  in  transit.  AIS  data  are  available  in  real  time  and 
document  the  actual  paths  taken  by  vessels  but  lack  the  cargo  details 
available  from  waterborne  commerce  data  though  CPT.  Thus,  the  arrival 
process  mining  techniques  discussed  in  this  report  serve  to  complement 
the  CPT  for  providing  USAGE  practitioners  highly  detailed  project-specific 
vessel  traffic  information.  Together,  CPT-  and  AlS-derived  information 
help  to  complete  the  waterfront  operational  picture. 

AIS  data  provide  a  practical  alternative  to  some  traditional  vessel 
observation  methods  as  spatial  and  temporal  information  is  captured  in  a 
single  technology  platform  and  avoids  some  problems  encountered  with 
traditional  observation  techniques  (Scully  and  Mitchell  2015).  Winkler  et 
al.  (2003),  Briggs  et  al.  (2004),  and  Maynord  (2007)  are  examples  of  field 
collection  efforts  that  would  benefit  from  use  of  AIS  data.  For  determining 
vessel  arrival  patterns  and  other  behaviors,  benefits  over  standard  practice 
arise  from  the  high  granularity,  standard  format,  automated  collection. 


ERDC/CHL  TR-17-2 


24 


lower  acquisition  cost,  and  centralized  aggregation  of  AIS  data. 

Operational  benefits  include  continued  collection  in  low  visibility, 
darkness,  and  situations  when  vessels  appear  at  differing  locations  with 
the  channel,  all  of  which  were  challenges  described  by  Maynord  (2007). 
Primarily,  the  decreased  cost  and  increased  scope  of  data  availability  that 
enable  nationwide  investigation  of  relevant  topics  are  demonstrated  here. 
AIS  technology  overcomes  most  environmental-related  limitations  arising 
from  adverse  weather  conditions  or  poor  visibility.  However,  AIS  data 
collection  is  limited  in  some  cases  by  line-of-sight  obstruction  from  vessel 
broadcast  and  shore-side  receiving  stations,  which  must  be  considered 
when  using  AIS  data  for  detailed  analysis. 

When  arrivals  are  visualized  as  shown  in  Figure  6,  periods  of  higher  (or 
lower)  use  within  a  single  port  become  obvious.  Differences  in  usage 
patterns  between  ports  also  become  obvious  when  these  distributions  are 
compared.  Arrival  frequency  distributions  may  provide  insight  that 
enables  port  managers  to  make  better  informed  decisions  in  support  of 
operations. 

The  automated  nature  of  AIS  data  broadcast  and  collection  allows  for 
automated  repetition  of  desired  analyses  according  to  project-specific  or 
programmatic  goals.  Data  for  an  entire  port  can  be  obtained  from  a  single 
source  by  employing  methods  described  in  this  report.  Consider  that  this 
type  of  tidal  analysis  could  be  performed  year  over  year  to  monitor  changes 
in  vessel  behavior.  Alternatively,  tidal  analysis  could  be  performed  after 
harbor  improvements  to  validate  assumptions  made  regarding  the  use  of 
available  water  levels,  which  would  inform  future  design  efforts. 

By  investigating  the  distribution  of  vessel  transits  with  regard  to  tidal 
elevation,  it  becomes  possible  to  attempt  grouping  and  categorizing  ports 
based  on  the  usage  patterns  of  incident  vessels.  For  instance.  Figure  9 
depicts  the  proportion  of  vessels  calling  above  each  port’s  respective  high- 
tide  threshold  versus  the  percentage  of  vessels  calling  below  that  port’s 
low-tide  threshold.  While  the  average  of  the  fraction  of  low-tide  traffic  is 
similar  across  regions,  Atlantic  ports  have  a  range  of  T25  values  triple  that 
of  Gulf  ports  and  nearly  double  that  of  Pacific  ports.  Similar  observations 
may  be  made  with  regard  to  the  range  of  T75  values.  Given  that  vessels 
arriving  anywhere  are  in  theory  free  to  call  at  any  tidal  stage,  the  cause  and 
meaning  of  these  observations  require  further  study  to  be  determined. 
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Figure  9.  Visualization  of  fraction  of  vessels  arriving  above  vs.  below  respective  tidal  elevation  thresholds 
by  region.  Much  more  variation  is  evident  at  Atlantic  and  Pacific  ports  than  at  Gulf  ports. 
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Figure  lo  highlights  the  possibility  of  using  the  tidal  dependence  metric  to 
monitor  changes  in  vessel  distribution  across  the  tidal  cycle.  Emergent 
trends  may  be  used  to  inform  decisions  related  to  maintenance  or 
expansion  investment  priorities.  For  instance,  Savannah,  GA,  shows  an 
increasing  trend  in  TD  values  over  this  3-year  sample.  This  is  too  brief  a 
record  from  which  to  draw  strong  conclusions,  but  Savannah  is  presently 
undergoing  harbor  deepening;  the  port  is  well  documented  as  being  visited 
more  frequently  by  vessels  that  are  tidally  constrained  due  to  large  drafts 
relative  to  available  depth.  The  decreasing  trends  in  TD  value  at  Portland 
and  Los  Angeles,  on  either  end  of  the  spectrum,  are  also  interesting. 
Further  conclusions  of  the  relevance  or  stability  of  these  trends  can  only  be 
made  by  expanding  the  temporal  scope  of  investigation. 
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Figure  10.  TD  values  by  port  and  year  visualized  for  year-over- 
year  stability. 
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Figure  ii  shows  the  solution  space  of  the  TD  metric  for  all  values  of  T25,  T50, 
and  T75.  This  figure  makes  it  clear  that  the  metric  is  highly  sensitive  to 
reductions  in  mid-tide  traffic,  especially  when  the  T50  fraction  approaches 
10%  of  the  total  traffic  population.  Below  10%  of  vessels  operating  at  mid 
tide,  the  TD  metric  is  increasingly  sensitive  to  imbalances  between  the  high- 
and  low-tide  traffic  fractions.  Within  the  full  range  of  possible  T50  values, 
the  raw  T25  and  T75  values  inform  cases  where  high  and  low  tide  are 
approximately  equal. 
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Figure  11.  Numerical  solution  for  possible  of  TD.  All  ports  investigated  for  this  work  fall  within  the  range  of 

0.3<T5o<0.6. 


Among  the  observed  ports,  the  mid-tide  portion  of  traffic  ranged  from 
0.46  at  Anacortes,  WA,  to  0.57  at  Boston,  MA.  Within  this  range  of  mid¬ 
tide  percentage,  the  maximum  expected  magnitude  of  TD  values  range 
from  +/- 1.17  to  +/-  0.75.  The  minimum  and  maximum  calculated  values,  - 
0.20  at  Long  Beach,  CA  (2013  and  2014),  and  0.55,  at  Portland,  ME 
(2012),  respectively,  fall  within  this  range.  The  TD  metric  as  currently 
formulated  is  limited  by  the  fact  that  the  variability  of  resultant  values  is 
perhaps  too  narrow  to  make  meaningful  comparisons  across  ports  with 
T50  values  near  0.5,  especially  when  T25  and  T50  values  are 
approximately  equal. 

Note  the  vessel  arrival  behavior  with  regard  to  the  value  of  the  TD  metric. 
The  maximum  individual  and  average  values  (0.55  and  0.39,  respectively) 
were  observed  at  Portland,  ME.  Portland  has  nearly  triple  the  volume  of 
arrivals  above  the  75th  percentile  water  elevation  (33%)  compared  to 
arrivals  below  the  25th  (13%)  on  average.  Conversely,  Los  Angeles,  CA, 
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demonstrated  the  lowest  average  TD  value  of  -o.i8  and  had  only  45% 
more  low-tide  calls  (29%)  compared  to  arrivals  at  high  tide  (20%).  Long 
Beach,  CA,  had  the  lowest  individual  TD  value  of  -0.20  with  traffic 
composition  similar  to  Los  Angeles.  Mid-tide  calls  were  nearly  evenly 
distributed  at  Portland  and  Los  Angeles,  with  53%  and  51%,  respectively, 
while  Long  Beach  had  50%  traffic  during  mid-tides. 

The  arrival  frequency  at  Los  Angeles  and  Portland,  shown  in  Figure  6, 
indicates  that  temporal  clustering  is  much  stronger  at  Los  Angeles  than  at 
Portland.  As  average  tidal  frequency  and  the  24  hr  day  are  unequal,  it  is 
expected  that  vessels  taking  advantage  of  high  tide  would  show  weak 
clustering  behavior  due  to  the  shifting  time  of  tidal  events  compared  to 
regular  daily  schedules.  Instead,  vessel  arrival  patterns  would  mimic  the 
occurrence  pattern  of  high-tide  events.  At  Los  Angeles,  where  call  clustering 
is  strongly  apparent  and  indicative  of  regular  daily  operations,  low-  and 
high-tide  arrivals  are  much  more  balanced.  This  is  to  be  expected,  again 
because  tidal  frequency  is  out  of  sync  with  a  24  hr  daily  schedule.  The 
interarrival  time  distribution,  shown  in  Figure  4,  indicates  that  vessel  calls 
are  spread  out  at  Portland  and  are  very  rapid  at  Los  Angeles  (which  had 
approximately  nine  times  as  much  traffic).  This  further  suggests  a 
relationship  between  traffic  volume,  temporal  clustering,  and  the  TD  value. 
Still,  the  preference  for  low  tide  at  Los  Angeles  is  curious  given  the  strong 
clustering  pattern. 

It  is  likely  that  the  variation  in  vessel  arrival  behavior  is  related  to  the 
vessel  types  and  the  drafts  at  which  those  vessels  arrive.  Figure  12  displays 
the  fractional  value  of  cargo  moved  through  Portland,  ME,  and  Los 
Angeles,  CA,  from  2012  through  2014  with  respect  to  draft,  as  reported  in 
Waterborne  Commerce  Statistics  Data  accessed  via  CPT.  The  volume  of 
throughput  is  substantially  larger  at  Los  Angeles  (60  million  tons, 
annually)  than  at  Portland,  (11  million  tons,  annually).  Los  Angeles  has  an 
authorized  entrance  channel  depth  of  53  ft  below  MLLW  but  logged  vessel 
arrivals  with  maximum  drafts  of  57  ft.  Portland,  with  an  authorized 
entrance  channel  depth  of  45  ft  below  MLLW,  logged  vessel  arrivals  with 
maximum  drafts  of  52  ft  during  the  study  period. 
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Figure  12.Vessel  traffic  at  Portland,  ME,  has  a  significant  value  of  cargo  moving  with  vessels  at  its  deepest  drafts. 
At  Los  Angeles,  CA,  the  value  of  cargo  transported  at  the  deepest  drafts  reported  in  the  harbor  is  marginal. 
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In  both  cases,  commodity  data  include  vessels  with  drafts  greater  than  the 
respective  authorized  channel  depths,  indicating  reliance  on  additional 
channel  depth,  very  likely  resulting  from  tidal  water  level  fluctuations.  The 
average  range  of  tidal  predictions  for  2012-2014  was  8.88  ft  in  Los 
Angeles  and  13.84  ft  in  Portland.  In  Los  Angeles,  additional  depth  from 
tide  accounts  for  16.8%  of  the  authorized  channel  depth.  In  Portland,  tidal 
depth  provides  an  additional  30.8%  beyond  authorized  depth. 

The  draft  of  shipments  in  excess  of  project  depth,  4  ft  at  Los  Angeles  and 
7  ft  at  Portland,  represents  7.5%  and  15.6%  of  authorized  depth, 
respectively.  The  draft  in  excess  of  project  depth  represented  45%  of  the 
observed  range  of  tidal  predictions  at  Los  Angeles  and  51%  at  Portland. 

Portland’s  cargo  is  predominantly  bulk  products  (99%  by  tonnage. 

Figure  13)  and  is  heavily  weighted  toward  its  deepest  drafting  vessels. 
Approximately  40%  of  cargo  value  moves  at  drafts  larger  than  the 
authorized  project  depth,  indicating  that  in  all  likelihood  it  transited  when 
water  levels  were  above  MLLW.  However,  less  than  1%  of  cargo  value  at 
Los  Angeles,  which  is  predominantly  containerized  products  (77%  by 
tonnage.  Figure  14),  moves  on  vessels  with  drafts  larger  than  the 
authorized  project  depth.  Based  on  AIS  data,  on  average,  33%  of  traffic 
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moved  above  the  75th  percentile  predicted  tidal  elevation  at  Portland.  At 
Los  Angeles,  only  20%  of  traffic  moved  above  the  75th  percentile  predicted 
tidal  elevation.  While  each  port  is  using  approximately  half  of  the  available 
tidal  depth,  Portland’s  increased  channel  depth  as  a  fraction  of  authorized 
channel  depth  is  approximately  twice  that  of  Los  Angeles.  The  benefit  of 
using  that  additional  depth  is  much  higher,  in  tonnage  and  value  terms  at 
Portland  than  at  Los  Angeles. 


The  preference  for  any  tide  stage  may  have  multiple  interpretations.  For 
example,  the  preference  for  higher  tide  elevations  at  Portland  may  reflect 
opportunism  on  the  part  of  vessel  operations  incented  to  land  the 
maximum  draft  possible  there.  Given  that  more  than  50%  of  traffic  in 
Portland  moves  during  mid-tide  elevations,  this  interpretation  seems  more 
likely  than  the  alternative  interpretation  that  bulk  cargo  vessel  operators 
might  prefer  to  avoid  tidal  currents  associated  with  mid-tide  elevations.  In 
Los  Angeles  temporal  clustering,  low  traffic  density  during  high  tide  and 
limited  value  of  cargo  moved  at  drafts  in  excess  of  authorized  project 
depth  may  mean  that  vessel  schedule  and  crane  productivity  dominates 
vessel  arrival  behavior  with  available  water  depth  being  generally  less 
important.  However,  detailed  analysis  of  the  motivating  factors  underlying 
these  general  observations  are  beyond  the  scope  of  this  investigation. 


Figure  IS.Portland,  ME,  cargo  is  predominantiy  buik  commodities.  On  average,  neariy  14%  of  cargo  moved  at 

52  ft  of  draft,  the  deepest  recorded  during  this  period. 
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Figure  14.  Los  Angeles,  CA,  cargo  is  predominantly  containerized.  Approximately  1%  of  cargo  on  average 
moved  at  the  deepest  recorded  draft  of  57  ft  during  this  period. 
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6  Conclusion 

This  investigation  has  identified  AIS  data  as  an  information  source  valuable 
for  collection  of  vessel-related  movement  information.  Provided  that  vessels 
of  interest  can  be  validated  through  authoritative  means,  AIS  provides  a  far- 
reaching  remote-sensing  platform  that  will  enable  practitioners  to 
efficiently  analyze  vessels  in  transit.  The  high  dimensionality  of  AIS  data 
provides  numerous  methods  for  filtering  data,  making  it  possible  to  analyze 
vessel  movements  precisely  at  an  unprecedented  breadth  of  scope. 

This  investigation  further  demonstrated  the  TD  metric  methodology  as  a 
simple  and  intuitive  measure  of  vessel  performance  within  USACE- 
maintained  navigation  entrance  channels  when  vessel  TD  imbalances  are 
present.  It  provided  a  relative  ranking  of  ports  within  the  USAGE 
navigation  portfolio  based  on  when  vessels  arrive  with  respect  to  predicted 
tidal  elevation.  When  extreme  (high  or  low)  tidal  preference  is  present  but 
balanced,  the  traffic  fractions  T25,T5o,T75,  are  more  informative. 

The  fusion  of  AIS  data  with  tidal  predictions  demonstrates  one  of  many 
potential  analyses  coupling  vessel  and  environmental  forcing  data.  AIS  is 
an  important  data  source  to  consider  when  investigating  other  navigation- 
related  topics,  especially  when  the  behavior  of  vessels  in  transit  is  relevant 
to  the  problem  under  consideration. 

Finally,  the  spatio-temporal  nature  and  high  granularity  of  AIS  data  make 
further  integration  of  vessel  transit  information  with  commerce  and 
channel  condition  data  a  worthwhile  area  for  future  research  and 
development.  In  the  case  of  Portland,  ME,  the  highest  observed  TD  value 
was  associated  with  a  substantial  value  of  port  cargo  being  moved  at  vessel 
drafts  in  excess  of  authorized  channel  depths.  In  Los  Angeles,  the  lowest 
observed  TD  value  was  associated  with  temporally  clustered  vessel 
arrivals.  While  both  observations  are  interesting  in  their  own  right,  it  was 
demonstrated  that  AIS  analysis  methods  serve  to  complement  existing 
tools,  such  as  CPT,  in  providing  USAGE  practitioners  with  richly  detailed 
information  that  helps  to  complete  the  picture  of  navigation  operations 
and  that  may  assist  waterway  managers  in  making  complex  channel  design 
and  maintenance  decisions. 
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Appendix  A:  Port  Reference  Lines 
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Vessel  Call  Frequency,  Boston,  MA,  2012 
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Vessel  Call  Frequency,  Boston,  MA,  2014 
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Boston,  MA  Vessel  Arrival  Water  Surface  Elevation 


Port  of  Interest:  Charleston,  SC 


Tide  Station  Number:  8665530 
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Vessel  Call  Frequency,  Charleston,  SC,  2012 
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Vessel  Call  Frequency,  Charleston,  SC,  2013 
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Vessel  Call  Frequency,  Charleston,  SC,  2014 
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Charleston,  SC  Vessel  Arrival  Water  Surface  Elevation 
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Astoria,  OR 
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Esri,  D^Ltarme,  GEBCO,  NOAA  NGIjC,  ard  ether  ocmiritjutQfB. 


Port  of  Interest:  Columbia  River,  OR 


Tide  Station  Number:  9439040 


Year 

Tide  Range  (Ft.) 

Z25(Ft.) 

Z75(Ft.) 

T25 

T5O 

T75 
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Arrivals 
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Vessel  Call  Frequency,  Columbia  River,  OR,  2012 
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Vessel  Call  Frequency,  Columbia  River,  OR,  2013 
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Vessel  Call  Frequency,  Columbia  River,  OR,  2014 
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Columbia  River,  OR  Vessel  Arrival  Water  Surface  Elevation 


Port  of  Interest:  Delaware  Bay 


Tide  Station  Number:  8537121 
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Tide  Range  (Ft.) 
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Vessel  Call  Frequency,  Delaware  Bay,  2013 
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Vessel  Call  Frequency,  Delaware  Bay,  2014 
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Freeport,  TX 

★ 


USCG  Freeport,  TX 

A 
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Esri,  D^Lcrme,  GEBCO,  NCAA  NGDC,  ard  ctJier  ooriributcfs. 


Port  of  Interest:  Freeport,  TX 


Tide  Station  Number:  8772447 
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Vessel  Call  Frequency,  Freeport,  TX,  2012 


-3  5  1  2  2  3  3  3  3  2  6  19121411  9 

-4  3  6  9  131012  5  0  2  4  6  7  1010  6 

-3  6  3  10  6  5  6  4  7[l^3  5  8  1110  7 

3  725512122227  10  84 
|-412567648p^698948 
6  353352332218995 


9  6  6 


11  5 

lI  ^ 
6  ^4  5 

6  6  2 

3  4  6 


4  7  2  5 
7  2  4  7 

5  6  0  8 
5  7  3  7 
7  11  9  4 
4  2  4  1 


2 

4- 

5- 
5- 
5- 
7 


f  M  M  M  ^  7  7  ^  ^  ^  7  ^  ^  f  ^ 


Hour (UTC) 


Vessel  Call  Frequency,  Freeport,  TX,  2013 
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Vessel  Call  Frequency,  Freeport,  TX,  2014 
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Freeport,  TX  Vessel  Arrival  Water  Surface  Elevation 
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Honolulu,  HI 

★ 

Honolulu,  HI 
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Esri,  D^jcrme,  GEBCO.  WOAA  l-JeDC,  ard  ether  cjatiiribetefs 


Port  of  Interest:  Honolulu,  HI 


Tide  Station  Number:  1612340 
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Esri,  Drfjcrme,  GEBCG.  WCAA  l-JCDC,  ard  ether  oofiiributars. 


Port  of  Interest:  Long  Beach,  CA 


Tide  Station  Number:  9410660 
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Vessel  Call  Frequency,  Los  Angeles,  CA,  2013 


-5  13  3 
11  5  10 
1414  4 
191411 
141416 
20  9  8 


4  5 
3  3 
3  6 
0  4 

5  5 


7 

7 

2 

5 

6 


4  6  8  2079  4317  2  2  1  3  2 

5  6  6  13  5136  27  6  3  3  6  2 

6 


65  64 


7  7  15 
7  9  8  151964 
6 


^416  5  2  3  3 
47  9  6  4  10  3 
3  5  8  19  36  3510  2  4  4  3 


9  11^3^19 
5  14[654625 
7  13204941 
3  4  252825 
5  4 182920 


1  3  3  8  6  9  7  1432g4110  4  4  2  2  1432342812 
f  ^  ^  ^  7  ^  ^  ^  ?  1,61,2  2,51,8  ? 


Hour (UTC) 


MON 
TUE 
WED 
THU 
FRI 
SAT 
SUN 

012345678  9  10111213141516171819  202122  23 


Vessel  Call  Frequency,  Los  Angeles,  CA,  2014 
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Hour (UTC) 


Vessel  Call  Frequency,  Los  Angeles,  CA,  2013 


MON 
TUE 
WED 
THU 
FRI 
SAT 
SUN 

012345678  9  10111213141516171819  20  2122  23 


-5  13  3 
11  5  10 
1414  4 
191411 
141416 
20  9  8 


4  5 
3  3 
3  6 
0  4 

5  5 


7 

7 

2 

5 

6 


4  6  8  2079  4317  2  2  1  3  2 

5  6  6  13  5136  27  6  3  3  6  2 

6 


65  64 


7  7  15 
7  9  8  151964 
6 


^416  5  2  3  3 
47  9  6  4  10  3 
3  5  8  19  36  3510  2  4  4  3 


9  11163  5619 
5  14[654625 
7  13204941 
3  4  252825 
5  4 182920 


1  3  3  8  6  9  7  1432g4110  4  4  2  2  1432342812 
f  ^  ^  ^  7  ^  ^  ^  ?  1,61,2  2,51,8  ? 


Hour  (GMT) 


MON 
TUE 
WED 
THU 
FRI 
SAT 
SUN 

012345678  9  10111213141516171819  202122  23 


Vessel  Call  Frequency,  Los  Angeles,  CA,  2014 
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Vessel  Call  Frequency,  Mobile,  AL,  2012 
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Port  of  Interest:  New  Haven,  CT 


Tide  Station  Number:  8465705 
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Vessel  Call  Frequency,  New  Haven,  CT,  2012 
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Vessel  Call  Frequency,  New  Haven,  CT,  2014 
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New  Haven,  CT  Vessel  Arrival  Water  Surface  Elevation 
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Port  of  Interest:  New  York  &  New  Jersey 


Tide  Station  Number:  8531680 


Year 

Tide  Range  (Ft.) 

Z25(Ft.) 

Z75(Ft.) 

T25 

T5O 

T75 

TD 

Arrivals 

2012 
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0.24 
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4817 

2013 
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0.26 

0.00 

4537 

2014 

7.88 

1.01 

4.14 

0.25 

0.49 

0.26 

0.01 

4471 

Tide  Cycle,  NQAA  Station  8531680,  January  1,  2012 
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Vessel  Call  Frequency,  Mew  York  &  New  Jersey,  2012 
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Vessel  Call  Frequency,  New  York  St  New  Jersey,  2013 
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Vessel  Call  Frequency,  New  York  St  New  Jersey,  2014 
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New  York  &  New  Jersey  Vessel  Arrival  Water  Surface  Elevation 


Port  of  Interest:  Norfolk,  VA 


Tide  Station  Number:  8638610 


Year 

Tide  Range  (Ft.) 

Z25(Ft.) 

Z75(Ft.) 

T25 
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T75 
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Vessel  Call  Frequency,  Norfolk,  VA,  2012 
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Vessel  Call  Frequency,  Norfolk,  VA,  2013 
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Vessel  Call  Frequency,  Norfolk,  VA,  2014 
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Norfolk,  VA  Vessel  Arrival  Water  Surface  Elevation 


Port  of  Interest:  Pascagoula,  MS 


Tide  Station  Number:  8741533 


Year 

Tide  Range  (Ft.) 

Z25(Ft.) 

Z75(Ft.) 
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Arrivals 
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Vessel  Call  Frequency,  Pascagoula,  MS,  2012 
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Vessel  Call  Frequency,  Pascagoula,  MS,  2013 
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Vessel  Call  Frequency,  Pascagoula,  MS,  2014 
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Pascagoula,  MS  Vessel  Arrival  Water  Surface  Elevation 


Port  Everglades 
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Port  of  Interest:  Port  Everglades,  FL 


Tide  Station  Number:  8723214 


Year 

Tide  Range  (Ft.) 

Z25(Ft.) 

Z75(Ft.) 

T25 

T5O 

T75 

TD 

Arrivals 

2012 

3.62 

0.44 

1.76 

0.24 

0.50 

0.25 

0.01 

3115 

2013 

3.57 

0.44 

1.76 

0.26 

0.49 

0.25 

-0.03 

3081 

2014 

3.80 

0.43 

1.76 

0.26 

0.49 

0.25 

-0.02 

3053 

Vessel  Call  Frequency,  Port  Everglades,  FL,  2012 
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Vessel  Call  Frequency,  Port  Everglades,  FL,  2013 
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Vessel  Call  Frequency,  Port  Everglades,  FL,  2014 
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Port  of  Interest:  Portland,  ME 


Tide  Station  Number:  8418150 


Year 

Tide  Range  (Ft.) 

Z25(Ft.) 

Z75(Ft.) 

T25 

T5O 

T75 

TD 

Arrivals 

2012 
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Tide  Cycle,  NOAA  Station  8418150,  January  1,  2012 
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Vessel  Call  Frequency,  Portland,  ME,  2012 
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Vessel  Call  Frequency,  Portland,  ME,  2014 
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Vessel  Call  Frequency,  San  Francisco,  CA,  2012 
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Tide  Cycle,  NOAA  Station  9755371,  January  1,  2012 
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Vessel  Call  Frequency,  Savannah,  GA,  2012 
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Vessel  Call  Frequency,  Savannah,  GA,  2013 
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Vessel  Call  Frequency,  Savannah,  GA,  2014 
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Vessel  Call  Frequency,  Seattle  &  Tacoma,  WA,  2012 
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Vessel  Call  Frequency,  Seattle  &  Tacoma,  WA,  2013 
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Vessel  Call  Frequency,  Seattle  St  Tacoma,  WA,  2014 
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Tide  Cycle,  NOAA  Station  8726384,  January  1,  2012 
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